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24 

25 ABSTRACT 

26 We have found that O- atoms are a primary product in the irradiation of water-ice with 0.8 

27 MeV protons. This observation has implications in understanding the chemical reactions 

28 that occur to produce molecular oxygen (O 2 ) in such laboratory ices, as well as ices found 

29 on the surfaces of Ganymede and Europa, and the ice particles present in Saturn’s rings. 

30 We estimate that in irradiated water-ice O 2 can be formed at a lower limit of 0.07% by 

31 number relative to water and is in agreement with observations of the icy Jovian 

32 satellites. 

33 

34 MAIN TEXT 

35 There have been several recent laboratory investigations of the formation of molecular 

36 oxygen from irradiated water ice (Orlando and Sieger et al., 2003; Petrik et al., 2005; 

37 Teolis et al., 2006). The formation of oxygen is of particular interest in the formation of 

38 O 2 on the surface of the icy Galilean satellites (Spencer et al., 1995; Spencer and Calvin, 

39 2002) and also an O 2 atmosphere associated with Saturn’s rings (Tokar et al., 2005; 

40 Johnson et al., 2006). In both cases, O 2 is postulated to be formed from the radiolytic 

41 and/or photolytic destruction of H 2 O molecules in the ices of these objects. 

42 

43 Despite ongoing efforts from various research groups, there is still a void in our 

44 understanding of the chemical processes that form O 2 molecules in irradiated water ice. 

45 The problem is not helped by the weak electronic and vibrational transitions of O 2 that 

46 make detecting solid-phase O 2 by optical spectroscopic methods difficult. Consequently, 
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the detection of 0 2 is commonly performed by quadrupole mass spectrometry (QMS) of 
molecules sputtered during the irradiation process, or later released in temperature 
programmed desorption (TPD) experiments. 

A good review of the field can be found in Johnson et al. (2003). Here we shall present a 
brief summary of the literature published since then and at the same time describe each 
0 2 production model. The focus will be on the chemical steps involved rather than the 
kinetics of each model. Details of our experiments will then be given. 

The papers by Johnson et al. (2003), Orlando and Sieger (2003) and Sieger et al. (1998) 
culminated in the publication of Johnson et al. (2006) with what shall be referred to here 

as the J(2006) model. This model proposes that 0 2 is produced from the decomposition 
of H 2 0 to H 2 and O-atoms. 


H 2 0 + energy H 2 + O 


( 1 ) 


The O-atoms are trapped in the form of a stable precursor, possibly an H 2 00 complex, 
before a second excitation produces H 2 and 0 2 . 


h 2 oo h 2 + o 2 


( 2 ) 


Alternatively a non-thermal O-atom from a secondary dissociation may form 0 2 . 
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70 H 2 0*0 + O H 2 O + O 2 (3) 

71 

72 An alternative model recently proposed by Petrik et al. (2006), and which shall be 

73 referred to here as the P(2006) model, proposes that the stable precursor is HO 2 and that 

74 more steps are required than in the J(2006) model. First, H 2 O is dissociated into H and 

75 OH. 

76 

77 H 2 0 + energy -> H+ OH (4) 

78 

79 Next the OH, which is formed within the sample, migrates to the surface of the ice. 

80 Multiple OH’s then react to form H 2 O 2 and subsequently HO 2 . 

81 

82 OH + OH H 2 0 2 (5) 

83 H 2 0 2 + OH “> H0 2 + H 2 0 (6) 

84 

85 The hydroperoxy radical (H0 2 ) is then dissociated by an energetic excitation to form 0 2 . 

86 

87 H0 2 + energy -> H + 0 2 (7) 

88 

89 Both models predict very specific but different atomic and/or molecular intermediates, 

90 yet in the experiments from which these models are constructed, the researchers never 

91 actually identified any of the chemical intermediates. 

92 
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93 The direct detection of radiolytically produced O-atom in an ice sample poses many 

94 experimental challenges. Perhaps the best evidence for the production of O-atoms in ice 

95 in the literature is found in the work of Matich et al. (1993) that detected Herzberg 

96 emission lines of 0 2 in UV-irradiated water ice. The authors reasoned that the 0 2 was 

97 formed from O-atom recombination, although other reactions not involving O-atoms 

98 were possible. Here we present new laboratory results on the observation of O-atom 

99 production in H 2 0 ice using the detection of isotopologues of ozone in irradiated H 2 0 + 

100 i8 02 thin-film samples. 

101 

102 The experimental details are similar to those described in Cooper et ah (2006). In brief, 

103 we prepared gaseous mixtures of H 2 “o + 18 0 2 (6:1) in a vacuum manifold. Millipore 

104 water was freeze-pump-thaw cycled multiple times to remove dissolved atmospheric 

105 gases. The ls 0 2 (Isotec; purity of > 97%) was used without further purification. Blank 

106 experiments on irradiated pure ls 0 2 produced l8 0., and no other detectable scrambled 

107 isotopes. The H 2 0 + ls 0 2 gaseous mixtures were then deposited onto an aluminum 

108 mirror, pre-cooled to 1 0 K by a closed-cycle helium refrigerator. The samples were then 

109 warmed to 80 K, at ~2 K/min. The 80 K irradiation temperature was chosen to 

HO approximate ice temperatures on the Galilean satellites. An increase in the vacuum- 

1 1 1 chamber's base pressure when the pre-irradiated ice was warmed to ~30 K indicated that 

1 1 2 some of the 0 2 sublimed out of the sample. This observation is consistent with previous 

1 1 3 work (Loeffler et al. 2006). Due to the sublimation of some of the 0 2 , the exact H 2 0/0 2 

1 1 4 ratio in the iiradiated ice is unknown. The samples were then irradiated at 80 K with 0.8 
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1 1 5 MeV protons generated by a Van de Graaff accelerator. IR spectra were measured using a 

1 16 Nicolet 6700 Nexus spectrometer at 4 cm -1 spectral resolution. 

117 

1 18 Figure 1 shows spectra of a H2 I6 0 + 18 Oa sample before and after irradiation. As the dose 

119 increased two broad but shallow absorption bands appeared at 980 and 990 cm' 5 

120 associated with the formation of the 888 and 886 isotopologues, where 8 represents an 

121 18 0-atom and 6 represents a 16 Oatom. The 888 was produced from residually trapped 

122 18 02 molecules that were nearest neighbors. The irradiation of pure 1S 02 did not yield any 

123 measurable ozone isotopes containing 16 0 and there was no measurable amount of CO 2 

124 atmospheric contaminant in our sample, so the 16 0-atom in the 886 must have been 

125 produced from H 2 O. This was the first ozone mixed isotopologue observed as the dose 

126 increased, and it is thought to be due to the first reaction step shown in Figure 2, i.e. the 

127 combination of a radiolytically produced 16 0-atom from water with 1S C >2 that was trapped 

128 in the ice. The 888 species was gradually destroyed (as shown in Figure 2) with 

129 increasing dose by the replacement of 18 0 by radiolytically-produced 16 0 to form 886. 

130 

131 We considered the possibility that the source of 16 0 in these experiments could be i6 OH 

132 formed in the radiolytie destruction of H 2 O. An ozone molecule containing a l6 0-atom 

133 could then be produced by the following reactions. 

134 

135 ,6 OH + 18 0 2 -> H 16 0 18 0 18 0 -> H + 16 0 18 0 l8 0 (8) 

136 
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137 However, calculations by Varandas (2002) and Yu and Varandas (2001) show that 

138 reaction 8 probably does not occur to any great extent since the energy of the H0 3 

1 39 intermediate is only marginally below that of OH + 0 2 on the ground-state H0 3 potential 

140 energy surface. These workers have shown that the O + H0 2 and H + 0 3 states are 

141 accessible in the gas-phase when the reacting 0 2 and OH are in highly excited vibrational 

142 states. These excited states are likely to be very rapidly quenched in the cold ices 

143 discussed in the present work. 

144 

145 Radiolytic dissociation from a stable H0 3 precursor (Reactions 9a-c) could also be 

146 possible. 

147 


H0 3 + energy -*> H + 0 3 

(9a) 

H0 3 + energy -> 0 + H0 2 

(9b) 

H0 3 + energy -> OH + 0 2 

(9c) 


151 

152 

153 However, from the potential energy surface (Varandas, 2002), reaction 9c would be 

154 expected to dominate due to the small difference in potential energy. Reaction 9a would 

155 be the least significant due to the largest difference in potential energy of reactants and 

156 products. In addition, the H0 3 abundance in irradiated H 2 0 + 0 2 ices at 80 K is small 

157 compared with the amount at 10 K (Cooper et al. 2006). 

158 

159 
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160 Once the 886 ozone species is produced, further reactions can occur. Other radiolytic 

161 !(I 0-atoms may add to either end of the 886 molecule and displace an O-atom from the 

162 opposite end. Alternatively, an ozone isotopologue may decompose to a molecular 

163 oxygen isotopologue and an O-atom (16 or 18). A 16 0-atom originating from radiolyzed 

164 water may then react with the 0 2 isotopologue to form either the 686 or 866 species 

165 which are observed at 1003 and 1022 cm' 1 respectively. The 866 can then form 686 or 

166 666, and the 686 may form 866 which can feedback into the former reaction channel. The 

167 666 isotopologue is observed at 1035 cm' 1 . The positions of these isotopologues are 

168 summarized in Table 1. 

169 

170 Assuming that oxygen atoms cannot add in between two other O-atoms in an ozone 

171 molecule, the 868 isotopologue would need to form from the reaction between an 886 or 

172 686 molecule and an 18 0-atom. The absence of the 868 (-1016 cm* 1 ) isotopologue in the 

173 sample indicates that 18 0 is irreversibly lost to the water lattice once it is displaced. This 

174 suggests that the amount of i6 0 produced radio lytic ally from H 2 0, must dominate the 

1 75 amount of 18 0 in the sample. 

176 

177 In our experiments on pure H 2 0, we see no infrared evidence of 16 0-atom production via 

1 78 the formation of molecular oxygen because the oxygen absorption band is very weak. 

179 Ozone is also not produced in any detectable quantity because the amount of 0 2 , the 

180 precursor needed for 0 3 production, is too small. However, when 18 0 2 is added to pure 

181 water, as described above, it acts as a trap for l6 0-atoms via the formation of ozone 

1 82 isotopologues that can be detected spectroscopically. 
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183 

184 Using a band strength of 1.4 x 10 17 cm molecule' 1 (Smith et al. 1985) for the v 3 feature of 

185 0 3 . and assuming that it is the same for all isotopologues, we have calculated that there 

186 are 0.14% 16 0-atoms by number relative to H 2 0 (Table 2) at the highest dose of 9.8 

1 87 eV716-amu molecule. In pure water ice, without the presence of 18 0 2 to trap I6 0-atoms as 

188 an ozone isotopologue, it could be argued that the highly reactive 16 0-atom may react 

189 with water or another water radiation fragment (such as H or OH) before ever 

190 encountering a second radiolytically produced ,6 Oatom. The lifetime of a single 16 0- 

191 atom surrounded by H 2 0 in an ice lattice is unknown, but it may be stabilized via 

192 forming a complex with H 2 0 (Khriachtchev et al., 2000). However, if the ozone 

193 isotopologues produced in the present work were dissociated and the I6 0~atoms were to 

194 reform exclusively as I6 0 2 , then there would be 0.07% i6 0 2 by number, relative to H 2 0. 

195 This small amount of 0 2 in pure water is far below the detection limits of our 

196 spectrometer, but using the O3 tracer molecule we can detect i6 0~atoms that could 

197 otherwise form 0 2 . 

198 

199 We note that these abundance estimates are calculated at the highest radiation dose of 9.8 

200 eV/16-amu molecule used in the present experiments. At this dose, l6 0-atom production 

201 does not appear to be at steady-state (Figure 3) i.e. the yield of O-atoms has not reached a 

202 point where the O-atom abundance doesn’t change with increasing dose. However, 

203 determining whether the O-atom production is at steady-state is impossible in these 

204 experiments because we cannot measure all of the O-atoms. We can only observe the O- 

205 atoms that are present as an ozone isotopologue. Previous experiments from our 
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206 laboratory (Cooper et al, 2008) have shown that in H 2 0 + l6 0 2 ices that 0 3 production 

207 reaches steady-state by -5 eV/16-amu molecule. 

208 

209 It is common in radiation chemistry to present production rates as a yield, G, the number 

210 of molecules produced per 100 eV of energy absorbed. We have calculated an effective 

211 Gh+( 16 0) by plotting the number of 16 0-atoms measured in ozone isotopologues as a 

212 function of radiation dose (Figure 3). The slope of Figure 3 is 0.013, a value equal to 

213 Gh+(* 6 0)- Again however, this value represents a lower limit due to the O-atoms that 

214 cannot be directly measured in this experiment and is only an effective G-value as we are 

215 actually not measuring 16 0-atoms directly, but indirectly in the form of ozone 

216 isotopologues. 

217 

218 Estimates of the 0 2 abundance on Ganymede range from 0.1 - 1.0% (Calvin et al. 1996) 

219 to 1.4 - 4.2% (Hand et al. 2006). While our percentage abundance is lower than these 

220 estimates, our value represents a lower limit estimate that is in good agreement with the 

221 planetary observations. In these laboratory experiments we cannot account for the l6 0- 

222 atoms that remain trapped in the ice and do not react, or the l6 0 2 that is formed but is 

223 present at a level below our detection limits because of the extremely weak fundamental 

224 vibration, and also the i6 0 2 formed that is sputtered or desorbed out of the sample. In 

225 addition, other trapping mechanisms such as clathrates (Hand et al. 2006) or simultaneous 

226 irradiation and deposition (Teolis et al. 2006), may increase the amount of 0 2 trapped in 

227 an icy satellite. 

228 
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229 In terms of modeling the radiation processes that produce 0 2 in pure H 2 0, we have found 

230 that O-atoms produced in these experiments, support the J(2006) model. While it does 

231 not by any means conclusively prove the model is correct, it does provide added support 

232 via the observation that O-atoms can be produced directly from H 2 0 molecules. The 

233 P(2006) model does not consider the production of such O-atoms in the formation of 0 2 

234 Further research is needed to produce a robust and accurate model for the formation of 0 2 

235 in irradiated, pure water ice, but we have shown here that O-atoms, produced directly 

236 from water may play a part. The actual mechanism will likely require a multi-instrument 

237 approach to detect all the intermediate species and reactions. In addition, further 

238 laboratory work needs to be performed to understand the trapping mechanism of 0 2 on 

239 icy satellites. 

240 

241 Nevertheless, we have shown that O-atoms are produced directly from H 2 0 molecules in 

242 an H 2 0 ice. If these O-atoms recombine to form 0 2 in similarly irradiated water ice, a 

243 lower abundance limit of 0.07% by number relative to water could result. This value 

244 agrees with the observed 0 2 abundance on Ganymede (Calvin et al, 1996; Hand et al, 

245 2006). 

246 
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295 Table 1: IR band positions (cm l ) for the v 3 vibration of ozone isotopologues in the gas 

296 phase (Schriver et al. 1 996) and in the present work in water ice. 



298 

299 

300 

301 Table 2. Percentage abundance of ozone isotopologues and the equivalent percentage 

302 abundance of l6 0-atoms in irradiated H 2 16 0 + 18 0 2 ice at 80 K after a dose of 9.8 eV/16- 

303 amu molecule. 

304 



307 
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308 Figure 1 : The absorption bands of ozone isotopologues produced in an H 2 0+ Chice 

309 sample irradiated with 0.8 MeV protons. 

310 

311 Figure 2: The chemical pathway for the formation of each ozone isotopologue. This 

312 model assumes that 16 0 only reacts with one end of an O 3 isotopologue and results in the 

3 1 3 loss of an O-atom from the opposite end. 

314 

315 Figure 3: The dose dependent 15 0-atom production (measured as an ozone isotopologue) 

316 at 80 K in an H 2 O + l8 02 ice sample irradiated with 0.8 MeV protons. 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 


Cooper O-atom paper 06-24-09-3 doc 



16 


331 

332 


334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 


Figure 1 


Wavelength (^m) 



Wavenumber (cm' 1 ) 


Cooper 0-atom paper 06-24-09-3 doc 




346 


Figure 2 


347 


348 


349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 


,8 o 18 o i8 o 

18 o Is o+ ,6 o-> 18 0 18 0 16 0 

' 8 0 ,fi 0 16 0<-> ,6 0 1 

!6 o 16 o 16 o 


8 0 i6 0 


Cooper O-atom paper 06-24-09-3.doc 



18 


365 

366 Figure 3 

367 


368 



Energy Absorbed (10 18 eV/ cm 2 ) 


Cooper O-atom paper 06-24-09-3.doc 




